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Abstract 



The use of superconductors in the field windings of 
large synchronous machines displaces iron as a major 
element in the production of the magnetic field and 
promises smaller, less costly machines for a given rating. 
The possibility of replacing the iron in a transformer 
bank with a rotating superconducting field is discussed. 

By placing the primary and secondary three phase windings 
in separate armature windings built in epoxy and fiber- 
glass concentric shells around the superconducting, free- 
spinning field and allowing the angular alignment of one 
of the armatures to be varied, a unique transformer is 
devised. Unique features of the synchronous transformer 
are: the ability to separately control real and reactive 

power flow, thus providing power factor correction while 
simultaneously acting as a transformer, and the ability to 
transform rated real power between two power systems that 
are out of time phase, with zero (or at a specified level 
within a range) reactive power flow. 
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Synchronous machines are commonplace in the modern world 
and the theory of their operation is well developed. The 
application of superconducting technology to the field wind- 
ings of synchronous machines is a new and growing area of 
research [1]. It is expected that superconducting machines 
with ratings of 1000 KVA and greater may have as much as an 
order of magnitude reduction in weight and size over a con- 
ventional machine of the same rating with only about, 0.1 
per cent of machine rating required to power the refrigeration 
unit of the superconductor [2]. 

Recent work at M.I.T. anticipates that it may be pos- 
sible to manufacture and ship synchronous generators with 
superconducting field windings with ratings as high as 
10,000 MVA and that construction and operation of these 
machines may be more economical than current conventional 
machines [ 3] . 

These economic advantages are brought about in part by 
the very high flux densities with minimal power dissipation 
achievable with superconductors, which makes it possible to 
eliminate the magnetic circuit except possibly as an external 
shield. Thus an expensive magnetic-iron rotor forging and 
stator teeth may be eliminated, which in turn reduces insula- 
tion requireme.n'f'S (or allows an increase in voltage rating) 
and at the same time increases the amount of space that may 
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be filled with armature conductors. Elimination of most of 
the iron from the machine reduces per unit reactances which 
may improve stability. 

References [2] and [4] derive expressions for electrical 
parameters of synchronous machines without magnetic shielding. 
References [5] and [6] extend these equations to machines 
with an exterior iron shield. 

In a superconducting alternator the superconducting 
field winding displaces iron as a major element in the pro- 
duction of the magnetic field. It might be asked if this 
can be done for other devices, specifically, is it feasible 
to replace the iron in a transformer bank with a rotating 
superconducting field? Such a device, hereafter called a 
synchronous transformer, would have two armature windings 
(one for each of two sets of three phase terminals) of the 
form contemplated for a superconducting alternator. It would 
also have a superconducting field winding, which is free 
spinning as does the field winding of a synchronous condenser. 

Electric energy is transferred from one armature to the 
other by the interaction of the two armature flux waves with 
the field winding flux wave: one interaction being similar 

to that of a motor, the other like a generator. 

As it will be shown, this proposed device has similari- 
ties with conventional transformers, but it has important 
differences too. For example, if the angular alignment of 
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one of the two armatures is made variable then the device 
has two "handles" on power flow: armature position and main 

field current. This contrasts with the availability of 
only one "handle" (tap changing) with conventional trans- 
formers. This thesis will show that these two control 
variables will enable the synchronous transformer to control 
separately both real and reactive power flow. Thus it can 
act as a power factor correction device while simultaneously 
acting as a transformer. Additionally, it can transform 
rated real power between two power system buses that are out 
of time phase with no reactive power flow or within a range 
of desired reactive power levels. 

Such a device could be constructed by embedding the 
armature windings in epoxy and fiberglass to a fiberglass 
(or other non-magnetic material) pipe for structural strength 
[7,6]. Two such armature windings would be required. The 
axis of the outer armature would be adjustable so that it 
could be rotated about the longitudinal axis of the rotor 
(and thus also about the longitudinal axis of the inner 
armature) by an external control, but it would not be free 
to spin. The returning flux path and shielding would be 
provided by a concentric iron shell external to the two 
armature shells. Figure 1 illustrates the physical configu- 
ration of the synchronous transformer. Figure 2 illustrates 
the electrical configuration. 
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Figure 1 

Physical Configuration of the Synchronous Transformer, 

End View 
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Figure 2 

Electrical Configuration 

Armature 2 is Shown Rotated by an Angle 4> from Armature 1 
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There is a growing interest in the use of supercon- 
ducting machinery for marine propulsion. In 1966 Sibley 
et al [8] estimated some potential advantages of super- 
conducting electrical machinery over conventional propulsion 
plants. In a later work Laine [9j considered the possibil- 
ity of using a superconducting transmission link in a 
destroyer- type ship. Matjasko [10] investigated a control 
system for a destroyer propelled by superconducting electric 
machinery. Greene [11] examined transmission systems 
including a cycloconverter scheme between a synchronous 
generator and an 18,000 SHP (13.5 MW) synchronous motor 
driving a propeller. 

While in its present configuration it would seem that 
the synchronous transformer would perhaps be more useful to 
the electric power industry than to marine propulsion, it 'is, 
however, similar in configuration to a class of dual arma- 
ture superconducting machines which may have great potential 
as an electric propulsion ship system. This class of AC 
machines may have the control characteristics of a DC machine 
system but with the additional advantage that a solid state 
switching circuit will replace the bulky cycloconverter. 

This class of dual armature machines will not be discussed 
further here as they are the proprietary information of the 
Cryogenic Engineering Laboratory at M.I.T. 
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Note : 


GLOSSARY OF TERMS 

Unless otherwise noted, the following subscripts 
will serve as modifiers to the symbols. 



Subscripts 



a 


phase a, armature 1 


b 


phase b, armature 1 


c 


phase c, armature 1 


a 


phase a, armature 2 


6 


i 

phase 6, armature 2 


Y 


phase y , armature 2 


d 


direct axis, armature 1 


q 


quadrature axis, armature 1 


6 


direct axis, armature 2 




quadrature axis, armature 2 


f 


field 


o 


zero sequence, armature 1 




zero sequence, armature 2 

% 


^1 


phase variables , armature 1 


^2 


pnase variables, armature 2 


^1 


rotor variables , armature 1 


^2 


rotor variables, armature 2 


r— 1 


armature 1 to field 


CM 


armature 2 to field 


^1 


equ: valent at the terminals , armature 1 


^2 


equivalent at the terminals , armature 2 



t 



Glossary of Terms 
( continued) 



current 
voltage 
winding flux 
inductance 

angular position of the rotor from the axis o 
phase a, armature 1, positive in direction of 
rotation 

angular position of the axis of phase a, 
armature 2 , from the axis of phase a, arma- 
ture 1, positive in the direction of rotation 

maximum value of mutual inductance between 
phase a, armature 1, and phase a, armature 2 
(i.e.,L at(J) = 0) 

maximum value of mutual inductance between 
phase a, armature 1, and the field 
(i .e. , L ^ at 0=0) 

maximum value of mutual inductance between 
phase a, armature 2, and the field 
(i.e. , at 0 = 1 ^) 

spatial angle from the rotor axis to the arma 
ture flux wave of armature 1, positive in the 
direction of rotation. 

spatial angle from the armature flux wave of 
armature 1 to the armature flux wave of 



armature 

rotation 


2, positive 


in the 


direc 


tion 


of 


terminal 


voltage 


of 


armature 


1 in 


per 


unit 


terminal 


voltage 


of 


armature 


2 in 


per 


unit 



per unit reactance 



air 



gap voltage 
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RESULTS OF A STATIC ANALYSIS 



For any winding of a synchronous machine the voltage at 
the terminals is equal to the product of the winding resis- 
tance and current through the winding plus the time rate of 
change of the winding flux. For any winding the flux is the 
summation of the self inductance times the current through 
the winding and the product of mutual inductances with 
other windings (which may depend on rotor angular position) 
and the currents through those windings. 

Thus : 



V 

m 



dX 



r 1 + o. 

mm dt 



m 



( 1 ) 



X 



m 



1 A 



i 

nm n 



(2) 



Losses in large synchronous machines are less than one 
per cent of machine rating and the assumption of a lossless 
model does not have much effect on the static analysis. 
Therefore one can assume that winding resistances are 
neglible yielding; 
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m 



dX 
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dt 



(3) 



by : 



The machine flux current relationships may be expressed 
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fixed 


geometry the following relations 


apply : 
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It is important to note that the geometry betv;een arma- 
tures is not fixed since armature 2 may be rotated and since 
phase axes are 120° apart. 



L 7^ L „ 7^ L 
aa a3 ay 



(9) 



It is assumed that all rotor to armature mutual 

inductances vary with the cosine of rotor position, with the 

maximum values of L , and L - when the rotor axis aligns with 

af af ^ 

the respective phase axes of armature 1 and armature 2. This 
is the same as neglecting higher order space harmonics in the 
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rotor air gap flux and assuming a non-salient rotor. 

It is further assumed that armature 2 to armature 1 mutual 
inductances vary as the cosine of the relative angular posi- 
tion of the respective phase axes with the maximum value of 
when the armature displacement angle, <p , equals zero. 

This value (maximum) will be denoted by M. This assump- 

tion is equivalent to neglecting higher order space harmonics 
in the air gap flux between armatures. Higher order space 
harmonics die away more rapidly than the fundamental.^ Since 
this is an air gap machine, disregarding higher order space 
harmonics in the air gaps will introduce only a small error 
[6]. Substitution into equation (4) yields equation (10). 



Equation (3) gives 
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V, 
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V, 




( 11 ) 



At this point it is convenient to transform equation 
(11) into the rotating coordinate system of the rotor. The 
subscripts d and q correspond to quantities of stator 1 , 
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rotating with the direct and quadrature axis of the field. 

The subscript, o, corresponds to zero sequence quantities of 
stator 1, that is in phase in all three windings of stator 1. 
The corresponding subscripts for stator 2 are 6, and 
The details of this transformation are given in Appendix A. 

The result of transforming to rotor coordinates yields 
equations (12) and (13) . 

By comparison of equations (12) and (13) the following 
relationships can be recognized: 



^d 


dt 


(A^) 


1 

E 

>* 


(14) 


V 




(A ) 


+ uA , 


(15) 


q 


dt 


q 


d 






il 


(A^) 


- o)A^ 


(16) 




II 


(A^) 


+ wA^ 


(17) 


^f 


II 


^f 




(18) 



It can be seen that zero sequence variables can be 
eliminated from equations (12) and (13) since they are com- 
pletely uncoupled from other variables. 

For ease in notation the following definitions are 
introduced : 



L„ — L — L . 
E a ab 



L = L - L i- 

e a aC 
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Thus equation (13) becomes: 
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Now that the basic relations governing the synchronous 
transformer have been developed, the operation of the trans- 
former between two buses will be looked at. Motor sign con- 
vention will be observed throughout, that is, real power 
input to the machine from either system (bus) is defined as 
positive and lagging current determines the positive direc- 
tion for reactive power. 



Since the flux-current relationships have been trans- 
formed into rotor d-q coordinates, the following definitions 
of armature equivalent flux vectors are introduced: 







(22) 




^t2 = ^^6" 


(23) 


Similarly, the 


equivalent armature voltage vectors 


are 


likewise defined: 








= •'''a' V 

T7 — VT7 ^ ^ 
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The angle 6^ is defined as the spatial angle from the 
rotor axis to the armature flux of armature 1, positive in 
the direction of rotation. The angle 6^ is defined as the 
spatial angle from the armature flux wave of armature 1 to 
the armature flux wave of armature 2, positive in the direc- 
tion of rotation. Figure 3 illustrates these sign conven- 
tions . 

Alternatively, since the voltage vectors are a ninety 
degree rotation of the armature flux wave vectors, these same 
definitions can be viewed as applying equally between the air 
gap voltage, terminal voltage of armature 1 and terminal 
voltage of armature 2. Figure 4 illustrates the sign conven- 
tions of angles 6^ and 5^ in terms of the voltage vectors. 

The following definitions follow directly: 



^d 


= “ V. 


. T sin6 
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cos6 


(27) 
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sin(6 + 6 ) 


(28) 
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t2 


a o 






= ^t2 


cos (6 + 6 ) 

3. O 


,(29) 



If the restriction to strict steady state operation is 
assumed, then time derivative terms become zero. Equations 
(14) through (17) become: 
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Figure 3 

Illustration of Sign Conventions 
for Spatial Angles 6^ and 6^ 
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Figure 4 



Illustration of Sign Conventions for Spatial Angles 
6^ and 6^ in Terms of Voltage Vectors 
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(33) 



Substituting of equations (21) through (29) into 
equations (30) through (33) yields: 

cos5^ = -cj(Lgig + M sin(j)i^ + M cos4>i^) (34) 

3 3 /3 

sin6^ = o)(Lgi^ + ^ cos<t>i2 - M sin(f>i^ + y (35) 

V.^ cos(6, + 6) = -cj(- 4 M sin(f)i, + 4 M cosdi + L i-) (36) 

0.0 z a z 

3 3 /T 

V. - sin(6 +6 ) = Lty{-^ M cos<f>i, + -;r M sin(f>i + L. ij + / rr L„i^ 

t2 ao 2 d2 q £0 /22f 

(37) 

Since the rotor is free spinning there will be no 
electrical-mechanical power transfer from the rotor. If we 
ignore all losses (friction/ windage, copper losses, etc.) 
then the net power of the synchronous transformer will be 
zero. 
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q q 
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(38) 



This is not to say that the device cannot absorb power 
from one bus under the lossless assumption, but rather 
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equation (38) tells us that real power into either 
armature must go out the other armature. 

Substitution of equations ( (14) through (17) with the 
steady state asumption d/dt = 0 and equation (21) into 
equation (38) yields: 
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o)L- i + 

1 q 






(39) 



Equations (34) through (37) and equation (39) are the 
governing equations for the synchronous transformer operating 
between two buses. Independent variables are i^, <P , 
and V^. Dependent variables are i^, i^, i^ , i^, 6^, and 5^. 
Thus we have five equations and six unknowns. The method of 
solution used is to assign values of 6^ for a sufficiently 
wide range and at each value of 6^ solve the resulting equa- 
tions which then have only five unknowns. However, before 
proceeding to the solutions of these equations it is con- 
venient to put them into tlie per unit system. Details of 
this are in Appendix B. 

The governing equations for the synchronous transformer 
operating between two buses in per unit are: 
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In order to use Gauss-Seidel computer aided iterative 
solution these equations were rewritten as follows; 
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V_ sin (6 +6 ) + X i, sind) - x i cosd) 
2 a o m d m q ^ 

X ~ 



(47) 



V» cos (6 +6 ) - X i, cosd) - x i sind) - x^ „if (48) 

2 ao md ^ mq ^ L2^ ' 



X , 



\ = - 



^L2^5 



X 



hi 



(49) 



Appendix (C) discussed the Gauss-Seidel iterative solu- 
tion. Reference 12 discusses this and other applicable 
methods in great detail. 

Reasonable values of per-unit inductances for a syn- 
chronous machine with a superconducting field winding were 
chosen [6, 13]. The assigned inductances are listed in 



Table 2. 
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PER-UNIT REACTANCES USED IN STATIC ANALYSIS 



Armature 1 Synchronous Reactance 



= 0.5 p. u 
= 0.5 p . u 
x^ = 0.4 p . u 
Xj^j^ = 1.0 p.u 
Xl 2 = 1*0 P*u 



Armature 2 Synchronous Reactance 



Mutual Reactance Between Armatures 



Armature 1 to Field Resistance 



Armature 2 to Field Reactance 



A computer program was prepared to solve simultaneous 
equations (45) through (49) by the Gauss-Seidel iterative 
method for a range of armature displacement angles, <p , from 
-90° to +90° and a range of spatial angles between armature 
voltage vectors, 6^, from 0° to 90° and for different levels 
of field current. In all cases reactances from Table 2 and 
one per unit terminal voltages were used. For each level of 
field current and each choice of angles and 6^ the computer 
calculated and printed values of direct and quadratic axis 
currents and voltages for both armatures , real and reactive 
power for both armatures , and the spatial angle between the 
air gap voltage vector and the voltage vector of armature 1. 
This produced a wealth of data from which Figures 5 through 
10 were extracted. 

It is important to note here that while the angles 6^ 
and 6^ are spatial angles relating the voltage vectors to the 
rotor, the anglf (}> is a physical angle by which the armatures 
are displaced.. Thus the time phase angle between the power 
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systems attached to the armatures is a combination of their 
physical and electrical angular difference, + <P- Since 
we are not interested in internal angles to characterize the 
device this power system time phase angle will be the sig- 
nificant angle used to characterize the device. 

Since we are primarily interested in the synchronous 
transformer to transfer power, the solutions yielding real 
power equal to one per unit were collated and graphed in 
Figures 5 through 7. Figures 8 through 10 were constructed 
from Figures 6 and 7. Figure 8 was obtained by cross plot- 
ting the armature displacement angle, (p , from Figure 6 
corresponding to the zero crossings of reactive power flow 
from Figure 7. Thus Figure 8 gives the setting points for 
armature displacement angle, <}), and field current. If, for a 
range of power system time phase differences in which one per 
unit real power can be transferred with no reactive power. 
Figures 9 and 10 yield similar results for reactive power 
flow equal to plus and minus 0.1 per unit. 



4> [in DEGREES] 
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So [IN degrees] 



Curves of armature displacement angle, <p , versus 
spatial angle between armature voltage vectors for 
one per unit real power transfer and various levels 
of field current, 




Curves of armature displacement angle, 4>, versus 
power system time phase angle difference, 6^ + 4*^ 
for one per unit real power transfer and various 
levels of field current, 









t 




Curves of reactive power, Q, versus power system 
time phase angle difference, 5^ + 4), for one per 
unit real power transfer and various levels of 
field current, 



- 2.0 




1.0 



0 : 

Ui 

0 . 






Figure 8 

Armature displacement angle, (p , and field current. 
If, versus power system time phase angle differ- 
ence, 6 q + 4> , for one per unit real power transfer 
and reactive power equal to zero. 



( 




-/ 00 - 



Figure 9 

Armature displacement angle, <t>, and field current. 
If, versus power system time phase angle differ- 
ence, 6q -1- 4), for one per unit real power transfer 
and reactive power equal to + 0.1 per unit. 



IN PER UNI T 



i 



36 







q: 






Figure 10 

Armature displacement angle, 4> / and field current, 
I^, versus power system time phase angle differ- 
ence, 6 q + 4), for one per unit real power transfer 
and reactive power equal to -0.1 per unit. 
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DISCUSSION OF RESULTS 



The important result of this thesis is that the 
synchronous transformer can be characterized by a set of 
curves in which the armature displacement angle and field 
current can be determined for a given power system time 
phase angle difference for operation at real power equal to 
one per unit and at various levels of reactive power. In an 
ordinary transformer there is always reactive power flow due 
to leakage reactance, but in the synchronous transformer 
settings may readily be chosen for unity pov;er factor power 
transfer. Additionally, the synchronous transformer can 
operate at unity power factor when transferring power between 
two systems that are out of time phase or it can simul- 
taneously act as a transformer and a synchronous condenser 
between two buses that may be in or out of time phase. Thus 
it is seen that a synchronous transformer in addition to 
being able to act as a transformer and as a synchronous 
condenser, can also act as a phase shifter. 
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CONCLUSIONS AND RECOMMENDATIONS 



Conclusions 

The results of this static analysis have been encourag- 
ing. It has been shown that the analytic model of equations 
(40) through (44) yield a set of curves in which the armature 
displacement angle and field current can be determined for a 
given power system time phase angle difference for operation 
at real power equal to one per unit and at various levels 
of reactive power. The synchronous transformer can transfer 
real power from one bus to another whether its time phase 
is behind, equal to, or ahead of that of the other bus. 

Recommendations 

It is recommended that further investigations of the 
synchronous transformer be pursued. Further efforts should 
include a dynamic analysis with faults and transients and 
an economic analysis to determine the best level of its 
usage within the power systems and its general economic 
feasibility . 
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The machine voltage current relationships (equations 
(10) and (11) may be expressed by equation (A-1) (page 43). 

Define a pair of reciprocal transforms relating phase 
and rotor corrdinates , where u represents a current or a 
voltage. These transforms are similar to those of Lewis [14] 



u 



u 



u 



= T 



u 



u. 



u 

I- c 



(A-2) 



u 



u, 



u 



= T 



-1 



u 



u 



^o 



(A-3) 





COS0 


cos(0 - 


2tt/3) cos ( 0 + 


27T/3) 




III 


-sin0 


-sin ( 0 


- 2tt/3) -sin(0 + 


27T/3) 


(A-4) 




_ 1/2 


1/2 


1/2 








COS0 


-sin0 


1 " 





T = 



/2 
/ 3 



cos (6 - 2 tt/ 3) -sin(6 - 2 tt/3) 1 

cos (6 + 2 tt/ 3) -sin (6 + 2 tt/3) 1 



(A-5) 



By partitioning equation (A-l) as shown we can reduce 
the notation as shown below and treat the partitions 



separately . 
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r 



»T3 

> 



XJ 

> 



I 

I 



1 
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H 

cu 

> 




A 

a 


A 

aa 


A ^ 
af 


^P2 


I! 


A 

aa 


A 

a 


^af 


1 

> 


1 


_ 


fa 


^f 



^P1 

^P2 



(A-6) 



Note that these subraatrices are not symmetric, i.e., 

A ^ . The off diagonal elements are transposes of their 

counterparts . 

The transformed version of equation (A-6) is: 



V 



R1 



V 



R2 



V, 



T;^:A T ^ 
dt a 



T;^A T ^ 
dt aa 



T;^A T ^ 
dt aa 



T;^ A T ^ 
dt a 



T;^ A 
dt af 



T;^ A 
dt af 



A. T ^ A. T ^ A. 

dt fa dt fa dt f 



■R1 



■R2 



(A-7) 



Carrying out the above matrix operations yields equa- 
tion (12) of the text. 

If the flux-current relationships (equation (10)) are 
partitioned in the same manner as equation (A-1) , they 
become : 




A 


A 


A 


a 


aa 




A 


A 


A 


aa 


a 




fa 


fa 


A 



af 

af 

f 



^P2 

if 



(A- 8) 



The transformed version of equation (A-8) is: 
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^R1 




T 


A T ^ 
a 


T A 

aa 


t“1 


T A ^ 
af 




^R2 


= 


T 


A t“^ 

aa 


T A 

a 




^ ^af 








A 


f 

fa 


A. t' 
fa 


-1 







(A-9) 

Carrying out this matrix operation yields equation 



(13) of thetext 
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There are various techniques for changing equations into 
per-unit equations, all of them related. The method that I 
shall employ parallels that of the notes for Power Systems 
Engineering II, course 6.552 at M.I.T., by Prof. G. L. 

Wilson . 

The following base quantities are specified: 



Axis 

d 

q 

6 

C 



f 



Base Voltage 

Y 

o 

V 



V 



N 
_C 

O N 



a 



V 



N 
_c 

o N 



a 



V 



o N. 



Base Current 



I 

I 

I 

I 



o 



o 



N 

o N 
N 

o N 



a 

a 

a 

a 



3 

2 



I 



N 

a 

o 



The subscripted N's in the above base quantities 
represent turns ratios. 

The equations to be put into per-unit are the governing 
equations of the synchronous transformer operating between 
two buses, equations (34) through (37) and (39) . These are 



listed below: 
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V, cos 6 
1 a 



3 3 

= - 0 ) (L„i_ + sin4> io + cose}) i^) 
Eq 2 o 2 5 



(B-1) 



V, sin 6 
1 a 



= o)(Lgi^ + |m cosci) - |m sin4> 

(B-2) 



cos (6+6 ) = -03 (- ^ sin4> i, + |m coscj) i + L i^-) (B-3) 

2 ao 2 <32 qe^ 



sin(6 +6 ) 
2 a o 



03 (|m cos<}> + |i4 sin<}> 



0 = 



03Liiq + 



(B-4) 

(B-5) 



However, we recall that equations (B-1) through (B-4) 
came from equations (14) through (17), which in steaciy state 
per-unit become: 



<''d>p.u. = 


0) 

Cjl) 

O 


(X ) 
q p. u . 


(B-6) 


<''q>p.u. = 


0) 


(X .) 

(3 p . u . 


(B-7) 


<''s)p.u. = 


0) 

0) 

o 


(X.) 

C P*u. 


(B-8) 


*'^5*p.u. 


Cl) 

CO 

O 


(Xp) 

6 p. u. 


(B-9) 



Thus we see that the quantities which we must fin<3 in 
per-unit are the fluxes. The flux equation is equation (21) , 
which is broken <3own into inc3ivic3ual flux equations below. 

* I “ ^6" I ^ /I hh 



(B-IO) 
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(A ,) 
d p . u , 



<V“o> 



(B-ll) 



0 ). 



^^d^p.u. V 



[ L„(i,) , I + 4 M cos4>(i.) I (N /N ) 

E d p.u. o 2 ^ 6 p.u. o a' a 



- ^ M sin4>(ii,) I (N /N ) + Lt (i,) (|- I ) (N /N,) ] 

2 ^ Cp*u. o a'a' /2 l f p.u. 2 ^ 

(B-12) 

3 3 

A = L„i + M sincf) M coscf) i- (B-13) 

cj ij q z 0 z ^ 



(A„) 



q'p.U. 



(B-14) 



0 ) ^ 

(A ) = [L„ (i ) I + 4 M sin4>(i^) I (N /N ) 

q p.u. Eqp.u. o 2 6p.u. oaa 






(B-15) 



Ag = I M coscj) + j M sintf) i + L^i + 



(A.) 



6^p.u. V N 

(— ) (— ) 
^to ^ 'n ^ 
o a 



(B-17) 






/I, 



N 



■*■ B (ic.) I (N /N ) + /„L„(i^) 



'e'-"5^p.u."o'^ a^^a' 2"^2 ^ p. u. '2 'n N 



f a 



(B-18) 



3 3 ' . . 

A^ = - -2 sin(j)i^ + -j M cos4>ig + 



(B-19) 
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(X^) 



c p.u. V N 

(— ) (— ) 
'■(0 ^ ^ 
o a 



(Xf) 

C P-u. 



(jJ T 

~ [ - ^ M sin4>(i,) I (N /N ) 
‘ 2 d'p.u. o a' a 



+ 4 M sin(})(i ) I (N /N ) 
2 ^ q'p.u.o a' a 



+ L (i^) I (N /N ) 2] 
e C p.u. o a' a 



Now identify the following quantities; 



(x ,) 

d p . u. 



CO L^I 
o E o 

V 



(x ) 
m p . u . 



^ CO MI N 
(1 _2 — 2.) (-2.) 
^2 V ^ ^ 

o a 



(x. ) 

6 p . u. 



CO L I N 
1 

o 



/ O £ O ' / a. 2 

'■ V ■' 'n ^ 



a 



^ --T (0 L, I N 

(yc \ =1/1 O 1 ° 

^ Ll^p.U. 2/2 V 

^ of 



- (0 L„I 

. _ 3 /3 o 2 o , a . 

^^L2^p.u. "2/2 V N^^ 

o ot f 



(B-20) 



(B-21) 



(B-22) 



(B-23) 



(B-24) 



(B-25) 



(B-26) 



. Now substituting equation (B-22) through (B-26) into the 
foregoing flux equations and dropping the parentheses and the 
subscript "p.u." since all quantities are now in per-unit 
yields ; 
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X, = x,i, + X if cosd) - X i- sin(t> + x^^i^ 
d dd mo m^ Ll f 


(B-27) 


X = Xji + X i. sind) + x i^ cosd) 
q dq m6 ^ m^ ^ 


(B-28) 


Xr = X i, cosd) + X i sind) + x-i. + x^ ..i^ 

6 md ^ mq ^ 6o L2f 


(B-29) 


Xr = ~x i, sind) + x i cosd) + x^-i- 
C md ^ mq ^ 6^ 


(B-30) 



Substituting these results into equations (B-6) through 
(B-9) and noting that for steady state synchronous operation 
w/co^ = 1.0 p.u. yields; 



V, = -(x,i + X i^, sind) + x i_ cosd)) 

d dq m6 ^ m^ 


(B-31) 


V = x,i, + x i. cosd) - x i^ sind) + x^^i, 

q dd mo ^ mC Llf 


(B-32) 


V. = x i, sind) - x i cosd) - 

6 md ^ mq ^ 6C 


(B-33) 


= x i , cosd) + X i sind) + x^ir- + _i, 

C md ^ mq ^ 66 L2f 


(B-34) 



Similarly, equation (B-5) can be put into per-unit as 
follows : 



0 = — L, (i ) I + — L-(i^) I (N /N ) 

(ij 1 q p.u. o 0 ) 2 5 p.u. o a a 

o o 


(B-35) 



Multiplication of both terms by the proper constant and 
assuming steady state operation yields: 



0 = T i + »i - 

Ll q L2 C 



(B-36) 



ii 

Equations (B-31) through (B-34) and (B-36) are the 
governing equations for the synchronous transformer in per- 
unit. 

The following definitions are introduced: 



V 



1 






u. 



(B-37) 



V 

2 



^ ^ • U 4 



(B-38) 



Substitution of the definitions of equations (26) through 
(29), (B-37), and (B-38), directly yields the per-uni,t equa- 

tions for the synchronous transformer operating between two 
buses . 



V,sin 6 = x,i + x i^ sind> + x i_ cos< 
1 a dq mo ^ mC 



(B-39) 



Vs (B-40) 



V 2 sin( 6 ^+ 6 ^) = -x^i^ sin4> + x^ig cos<|) + x^i^ (B-41) 



V_cos (6 +6 ) = X i, coS(j) + X i sin4> + x^^i. + x^„i_ (B-42) 

2 ao md ^ mq ^ 66 L2f 



= ^Ll^q + ^L2^C 



0 



(B-43) 
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GAUSS-SEIDEL METHOD 



The Gauss-Seidel method is a useful technique for 
solving the following system of simultaneous linear equations 

a,,x, + a,-X- + ... + a, X = a, 

11 1 12 2 In n l,n+l 

a^.x, + a„„x« + . . . + a„ X = a„ ,, 

21 1 22 2 2n n 2,n+l 

• • 

• • 

• # 

a ,x, + a -x„ + ... + a X = a , 

nl 1 n2 2 nn n n,n+l (C-1 

In order to start the algorithm an initial solution 

vector must be assumed and the interative algorithm is: 



a . , , n a . . 

X. = ^ ^ X. 

^ ^ii j = i "ii 

i¥± 



(C-2) 



In the iterative algorithm, (C-2) , the most recently 
available values of x^ are always substituted in on the 
right-hand side. 

The convergence criterion is: 

|x^(k+l^^ iteration) - x^ iteration) | < e (C-3) 

r 1, 2, ..., n 



For computer solution an upper limit on the number of 
iterations should be specified since convergence may not 



occur . 
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